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Abstract Overlap fermions preserve a remnant of chiral symmetry on the lattice. 
They are a powerful tool to investigate the topological structure of the vacuum of 
Yang-Mills theory and full QCD. Recent results concerning the localization of topo- 
logical charge and the localization and local chirality of the overlap eigenmodes are 
reported. The charge distribution is radically different, if a spectral cut-off for the 
Dirac eigenmodes is applied. The density q{x) is changing from the scale-a charge 
density (with full lattice resolution) to the ultraviolet filtered charge density. The 
scale-a density, computed on the Linux cluster of LRZ, has a singular, sign-coherent 
global structure of co-dimension 1 first described by the Kentucky group. We stress, 
however, the cluster properties of the UV filtered topological density resembling the 
instanton picture. The spectral cut-off can be mapped to a bosonic smearing pro- 
cedure. The UV filtered field strength reveals a high degree of (anti)selfduality at 
"hot spots" of the action. The fermionic eigenmodes show a high degree of local 
chirality. The lowest modes are seen to be localized in low-dimensional space-time 
regions. 
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1 Introduction: Overlap Fermions and Topological Charge 

Quantum chromodynamics (QCD) is the theory of strong interactions. It is formu- 
lated in terms of quarks and gluons. The task is twofold: to describe the composite 
structure and the high-energy interactions of strongly interacting hadrons, in both 
cases taking the substructure in terms of quarks and gluons into account. One prob- 
lem for theorists and experimentalists is that quarks are permanently confined mside 
hadrons. Apart from details, the spectrum and symmetries of hadrons are dictated by 
an approximate chiral symmetry and its spontaneous breaking by the interaction via 
gluons. These properties are unique for QCD as part of the standard model and re- 
sult from the vacuum fluctuations of gluons. At a temperature of about 160 MeV the 
hadronic world experiences a transition to a quark-gluon plasma phase with rather 
unusual properties. Then, the vacuum structure has changed. 

Simulations on a space-time lattice are the only ab initio approach to these phe- 
nomena. This approach has the virtue that also structural information on the vacuum 
fluctuations is accessible. On the other hand, there are models attempting to give a 
qualitative understanding. Some of them, like the instanton model, were partially 
successful. The instanton model became challenged more recently. A study of how 
the topological charge density is distributed in space-time and how the gluon field 
localizes the quarks (in analogy to the Anderson effect) is crucial to understand the 
microscopic mechanisms. 

Overlap fermions [1,2] possess an exact chiral symmetry on the lattice [3] and 
realize the Atiyah-Singer index theorem at a finite lattice spacing a [4] . This is possi- 
ble because they allow a clear distinction between chiral zero modes and non-chiral 
non-zero modes. Depending on their chirality, counting of n+ or n_ zero modes de- 
termines the topological charge of the gauge field as 2 = n_ — n+ . Furthermore, they 
give rise to a local definition of the topological charge density [5]. Altogether, this 
makes overlap fermions an attractive tool for investigating the chiral and topological 
structure of the QCD vacuum. 

The topological structure attracts attention because of the old hope that local ex- 
citations contributing to the winding number might not only reaUze the breaking 
of Ua(1) symmetry but simultaneously provide a mechanism for confinement and 
chiral symmetry breaking. The instanton liquid model (ILM) does not fulfill this 
expectation: it can account for chiral synraietry breaking but fails to give an expla- 
nation of confinement. 

When the QCDSF collaboration had started to analyze ensembles of lattice con- 
figurations with overlap fermions for quenched QCD [6, 7] and later for QCD with 
dynamical quarks [8] by diagonalizing the massless overlap Dirac operator and to 
store an - until then unprecedent - number of eigenmodes (#"(150) per configura- 
tion), the way was open for a serious investigation of the structure of topological 
charge [9-11 J. Not long before we started this investigation, the instanton model 
had been challenged [12-14] by the observation that the topological charge, rather 
than appearing in Ad clusters, possesses a global, sign-coherent ?>d membrane-like 
structure [15-17]. We shall discuss here what remains from the instanton picture. 

The overlap Dirac operator D has to fulfill the Ginsparg- Wilson [18] equation. 
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Y5D-^+D-^Y5=a2RY5, (1) 

with a local operator^. This is what maximally can remain of chiral symmetry on 
the lattice. A possible solution - for any input Dirac operator, i.e. for the Wilson- 
Dirac operator Dw as well - is the following zero mass overlap Dirac operator 

£)(m = 0) = ^|l+ .r:^ I =^(l + sgn(£);y)), (2) 



with Dw = M — ^ where M is the Wilson hopping term and ^ a negative mass term 
to be optimized. This operator is diagonalized using a variant of the Amoldi algo- 
rithm. The operator can be improved by projecting the Ginsparg-Wilson circle on 
the imaginary axis, A Aimp. The topological density can be defined (with maximal 
resolution a) as follows: 




q(x) = -tr 75 (l - ^£)(m = 0;x,x)) 



(3) 



Using the spectral representation of (3) in terms of the eigenmodes V'a(^) with 
eigenvalue A, an UV smoothed form of the density can be defined by filtering, 

<lx.M = - I fl-y)l(rAW'75rIW), (4) 

summed over color c and with Asm as an UV cut-off. Our intention was to study the 
topological charge density at all scales, in particular the spectral representation in 
terms of overlap modes defining qx^^{x). Also, we had the opportunity to evaluate 
the topological charge density without filtering (i.e. at the resolution scale a) on the 
Linux cluster of the LRZ in Munich. We call this the scale-a topological density 
q{x) in contrast to qx^^{x) defined with UV smoothing. (For details see Ref. [11].) 



2 Topological Density with Different Resolution 

Figure 1 shows the topological density calculated with four different smoothing 

scales Asm given in units of 1/a. The presentation is in the form of isosurfaces at a 
fixed value of \q{x) \ = q^ut- With increasing Asm> the number and shape of the visible 
clusters changes. The given configuration has Q=l. 

Therefore exactly one zero mode exists. In Fig. 2 we show two extreme cases, 
the highly localized contribution of the single zero mode alone (left) and the scale-a 
topological density (right). For the same value of ^cut the isosurfaces of the latter 
fill the whole volume, separating essentially a single positive from a single negative 
cluster. The net charge is naively expected to reside in an unpaired instanton (located 
at the zero mode). Actually, however, the positive and negative charge (differing by 
one) is globally distributed over the two clusters. 
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Fig. 1 Distribution of the topological charge density in a given time slice of a 1 1-* x 24 lattice at 
j8 = 8.10: the spectral cut-off a As„ =0.14, 0.28, 0.42, 0.56 (from left to right). The isosurfaces are 
shown for q{x) = ±0.0005. Colors red/green denote the sign of topological charge. Fig. from [9]. 




Fig. 2 The same plot as in Fig. 1 for the zero mode contribution (left) and the unfiltered density 
(right). Fig. from [9]. 

3 Cluster Analysis 

More quantitative information can be obtained from a cluster analysis. The value 
qcm, that has characterized the isosurfaces in Figs. 1 and 2, is now treated as a run- 
ning parameter. The region characterized by \q{x)\ > qcm (i.e. the interior of the 
isosurfaces) consists of some number of mutually disconnected clusters. We have 
classified them according to their size and other properties, for example their fractal 
dimension. Results of the cluster analysis are presented in Fig. 3. In Fig. 3a we show 
the number of clusters as function of ^cut/?max for three lattices characterized by dif- 
ferent coarseness. Except for the coarsest lattice with j3 = 8.1 (a = 0.142 fm), the be- 
havior is similar. Around g'cut/^max = 0.2 . . . 0.25 the multiplicity of clusters reaches 
a maximum with a cluster density of ^(75 fm"'*) for )3 = 8.45 (a = 0.105 fm). 
The density is even higher for the finer lattice with j3 = 8.60 (a = 0.096 fm). In the 
limit ^cut/^max the number of clusters reduces to the two oppositely charged 
global clusters. Figure 3b shows the volume fraction of the biggest cluster rela- 
tive to all clusters. Except for the coarsest (j3 — 8.1) lattice, above ^cut/^max = 0.2 
where the number of clusters is still growing towards the maximum, this fraction 
is small (< 10 %). It changes rapidly below q'cut/^max = 0.2 and approaches 50 % 
for ^cut/?max < 0.1 and remains so in the limit when only the two biggest clus- 
ters remain. The distance between the two largest clusters is shown in Fig. 3c. The 
"distance" is defined as the maximum taken over all points in one cluster of the min- 
imal distance to any point of the other cluster It approaches 2a for ^cut/^max < 0.1, 
meaning that the two clusters are closely entangling each other There are no points 
deep in the interior, such that the clusters must be considered to be of lower dimen- 
sion than 4d. Figure 3d shows the connectivity (explained in [11]) of the biggest 



Topological Structure of the QCD Vacuum Revealed by Overlap Fermions 



5 




Fig. 3 Cluster analysis of the unfiltered topological density: the (jcut-dependence is shown (a) of 
the number of separate clusters, (b) of the volume of the largest cluster relative to all clusters, (c) of 
the distance between the two largest clusters in lattice units and (d) of the connectivity describing 
cluster percolation [11]. The data is plotted for the 12^ x 24 lattice configurations at j8 = 8.10 (red), 
for the 16"' x 32 lattice configurations at /5 = 8.45 (green) and for the 16"' x 32 lattice configurations 
at J3 = 8.60 (blue). Fig. from [11]. 



cluster, describing the ability to span the whole lattice. The connectivity signals the 
onset of percolation which begins below ^cut/?max = 0.2, too. All this confirms the 
description of "melting instantons" given earlier by the Kentucky group [12-14]. 

The meaning of the cluster analysis is geometrically visualized in Fig. 4, again 
contrasting the scale-a density with the filtered density. Both types of clusters are 
presented as a function of <7cut- The upper row shows, say for q'cut/^max ~ 0.2, that 
more and more irregular, spiky clusters of the unfiltered density are popping up be- 
fore percolation sets in. The lower row demonstrates that the number of clusters of 
the filtered density remains small and practically independent of ^cut- The density is 
^(1 fm""*), that means comparable with what has been estimated for the instanton 
density. Of course, with lowering ^cut, the clusters shghtly grow before they eventu- 
ally merge. Coalescence happens finally below ^cut/^max = 0.1. We see that the UV 
smoothed density exhibits cluster properties similar to the instanton model [19] or 
what cooling studies on the lattice have shown earlier [20,21]. 

4 Fractal Dimensions 

The visualization in three dimensions does not immediately expose the (fractal) di- 
mension of the clusters of the scale-a topological density. We have measured the 





Fig. 4 Isosurfaces of topological charge density with \c/{x)\/q„t^x = 0.1, 0.2, 0.3, 0.4, in one times- 
lice of a single 16"* x 32 configuration at j3 = 8.45. The upper pictures are based on the scale-a 
density, the lower pictures on the eigenmode-truncated density with a Asm = 0.076. Color encodes 
the sign of the charge enclosed by the isosurface. Fig. from [10]. 



(fractal) dimension of these clusters by a random walk method. While the cut- 
off qcm defines the extension of the clusters as discussed above, a random walk 
is arranged inside the biggest cluster that is unable to penetrate the (fractal) bor- 
der. A rough estimate of the dimensionality of the cluster at the respective level 
of ^cut can be inferred from the return probability to the starting point of the ran- 
dom walk. In d* dimensions it behaves as a function of the number of steps like 
gure 5 shows the observed power like behavior for differ- 
ent ^cut/^max- The dimensions d* extracted from this study are shown in the legend. 
This analysis reveals a continuous change of the fractal dimension from d* < I 
(characterizing the irregular spikes) to t/* « 3 (characterizing the 3D membranes) 
with a lowering cut-off ^cut/?max — 0.4 . . .0.0. 
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Fig. 5 The return probability of random walks in the biggest cluster of the scale-a topological 
density as a function of the step number, and the fractal dimensions infeiTed from that for various 
lower cut-off values q^^t- Fig- from [11]. 
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5 Smearing vs. Filtering 

The similarity of the filtered topological density to clusters obtained earlier by 
smearing methods can be studied more in detail by direct comparison. First let us 
recall that the overlap operator is not ultralocal. Even the scale-a topological den- 
sity given by overlap fermions represents some inherent non-locality of the overlap 
operator [22]. This is seen by considering the two-point function of the topologi- 
cal density. Theoretically, one expects that it is negative for all non-zero distances, 
with an infinite, 5-like contact term at vanishing distances. The two-point function 
actually measured is shown as the green curve in Fig. 6a. One needs to go to larger 
/3 (a finer lattice) in order to see the positive core slowly shrinking and the mini- 
mum becoming deeper. To compare with smearing, the alternative measurement of 




X [!m] X [!m] 

(a) (b) 

Fig. 6 Left: the two-point function of the fermionic topological density at scale a (green) compared 
with the bosonic definition after 5 steps of stout-link smearing (blue). Right: the two-point function 
of the fermionic topological density with an UV cut-off at Asm = 634 MeV (red) compared with 
the bosonic definition after 52 steps of stout-link smearing (blue) when the correlator is fitted best. 
The best point-by-point matching of the density q^^^ with Asm = 634 MeV is achieved with 45 
steps. For this case of less smearing the correlator is steeper. Fig. from [23]. 



the topological density (and of the two-point function) uses an improved bosonic 
definition of the topological charge density, essentially 

16 

in terms of an highly improved electric and magnetic field strength [24], that in- 
cludes plaquettes and bigger loops. To apply it to a generic lattice configuration 
requires a few iterations of smearing (actually, of stout link smearing with respect 
to an overimproved action [25]) before integer-valued topological charges Q are 
measured with good precision. Figure 6a shows that the overlap definition and the 
bosonic definition of the topological density agree best to reproduce a similar two- 
point function at 5 smearing steps. Figure 6b shows the gluonic two-point function 
after 45 and 52 smearing steps compared with the overlap definition of the two-point 
function with Asm = 634 MeV. At this level of smearing, the negativity is already 
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washed out and extended clusters are visible. A point-by-point matching of the topo- 
logical density and a matching of the two-point function is achieved with almost the 
same number of smearing iterations. 

There are two parameters, ^niax and the size pinst (historically the instanton ra- 
dius) suitable to characterize each cluster. Now the size parameter pinst is estimated 
by the curvature of q{x) close to the maxima xq where q{xQ) — ^max- In the form 
of a scatter plot in the plane spanned by these two parameters one can describe and 
compare the cluster composition of the gauge field. In Fig. 7 this is shown after 5 
and 40 stout-link smearing steps. After 40 smearing steps the cluster multiplicity is 
essentially reduced, while an instanton-like relation (the curves drawn in Fig. 7) 



^max 



rinst 



(6) 



is enforced within some tolerance. At the same time the clusters are far from being 
ideal 0(4) symmetric instanton solutions. The bosonic and the fermionic view of 




p [lattice units] 




p [lattice units] 



(a) (b) 
Fig. 7 The cluster content after 5 and 40 smearing steps in the (/max-Pinst plane. Fig. from [23]. 



the cluster structure is shown in the left and right panel of Fig. 8. They show a 
remarkable similarity. 

6 Selfduality 

This leads to the question to what extent the UV smoothed fieldstrength tensor is lo- 
cally selfdual or antiselfdual, or even a semiclassical (instanton) solution. Gattringer 
was the first to ask this question in Ref. [26]. We have applied his technique using 
the eigenmodes of the overlap Dirac operator. Truncating the sum over eigenvalues, 
an UV filtered form of the field strength F^y can be obtained that allows to evaluate 
the respective infrared (IR), UV smoothed topological charge density, 

qiR{x) oc Tr {F^y (x) F^y (x)) , (7) 

and action density 
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(a) (b) 



Fig. 8 The fermionic topological charge density of a 2 = configuration with = 634 MeV 
(left) compared with the bosonic one after 48 sweeps of stout-link smearing (right). Negative: 
blue/green, positive: red/yellow. Fig. from [23]. 



SIR (x) oc Tr [F^v (x) F^v {x)) ■ (8) 
Analogously to the local chiraUty of the non-zero modes one can define the ratio 

siR[xj + qm[x) 

which can be converted to 

/;(x) = ■^arctan(^V7(];o) -1 e [-1,+!]. (10) 

Regions with R{x) ^ —I or R{x) « +1 are characterized by the field strength be- 
ing locally selfdual (antiselfdual). Figure 9 shows histograms with respect to R for 
the whole lattice and for subsets of lattice points selected by ranking according to 
the action density. With higher action density (say, for < 30% of the lattice points 




Fig. 9 Normalized histogram with respect to the UV smoothed (anti)selfduality R for different 
subsamples of lattice points (1%, 5%, 10% etc.) ordered with respect to the UV smoothed action. 
20 overlap modes have been included in the smoothing. Fig. from [11]. 
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forming "hot spots" of filtered action Sjr) the parameter /J is distributed close to ±1 
(i.e. \R\ > 0.7). Clusters with respect to \R\ > 0.98 are overlapping with clusters of 
a suitably filtered topological density. 

7 Localization and Local Chirality of Overlap Eigenmodes 



A Im lA,i,„pl = 50 MeV 
- p = 8.10, 12^24 
-H-- p = 8.45, 12^24 
-A- p = 8.45, 16^32 
p = 8.45, 24^48 
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Fig. 10 IPR: the dependence on A (left), the a-dependence (right upper), and the V-dependence 
(right lower). Fig. from [9]. 



We have also analyzed the localization behavior of individual overlap eigen- 
modes. The interest in the localization properties of Dirac and scalar eigenmodes 
has been discussed by de Forcrand [27]. In the left panel of Fig. 10 we show the 
average inverse participation ratio (IPR) typical for zero modes and for bins of 
eigenvalues in dependence on coarseness and lattice size. The localization grows 
approaching the continuum limit and with increasing physical volume. The right 
panel shows, separately for zero modes and the first bin ( | A | < 50 MeV), the change 
of the IPR with a and with V. From fits we have concluded that zero modes are typi- 
cally 2-dimensionally extended, while the lowest non-zero modes are 3-dimensional 
objects. 

Similar to (9) the local chirality X {x) of the non-zero modes can be defined. This 
quantity is highly correlated with the UV smoothed topological density. For Asm = 
200 MeV the correlation function between X{x) and qx^^^iy) is shown in Fig. 11. 
The correlator is positive for the lowest 120 eigenmodes within a range of distance 
R < ^(1 fm). 
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R 

Fig. 11 The correlation function between the local chirality X{x) of selected eigenmodes and the 
filtered topological density q;^ (y) in a configuration with topological charge 2 = 0. Fig. from [11]. 



8 Technical Details 

To compute the low lying eigenvalues of the overlap Dirac operator (2) we use an 
adaption of the implicitly restarted Arnoldi method which is used in the ARPACK 
[28] library. The advantage of this method is that the action of a matrix on a vector 
can be computed freely without the need to express the input matrix explicitly. 

The main computational challenge using the overlap Dirac operator as an in- 
put matrix is the computation of the sign function sgn{Dw) in (2). We project out 
the lowest ^(50) eigenvalues of Dw and treat them exactly. The rest acting in the 
orthogonal subspace is approximated using minmax polynomials [29]. The applica- 
tion of the Wilson-Dirac operator Dw on a vector, the main computational kernel of 
almost every lattice QCD code, is implemented in assembler and highly optimized 
using SSE2 instructions. The data is aligned on full cache lines and stored in a spe- 
cial way to make optimal use of the 128-bit XMM registers. Prefetch instructions 
are issued to move the data to the cache accurately timed. The SSE2 version reaches 
approx. 1/3 of peak performance on a 16^ x 32 lattice on the Linux cluster of the 
LRZ. 

The computation of the scale-a topological density q{x) requires to evaluate the 
local trace over Dirac and color indices of the overlap Dirac operator with the 75 
matrix included according to (3). This can be trivially parallelized, but is numeri- 
cally extremely expensive, since the overlap operator has to be applied on 12 Viat 
unit vectors. To compute q{x) on one single configuration on a 16^ x 32 lattice 
approx. 1.3-10^ CPU hours for totally 1 .6 million applications of the overlap Dirac 
operator were required. 
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9 Conclusions 



We have confirmed that the scale-a topological density is highly singular and the 
topological charge of either sign is globally filling the two extended 3d percolat- 
ing structures. A similar low-dimensionality is found for chiral zero-modes (two- 
dimensional) and the lowest non-chiral non-zero modes (three-dimensional). The 
microscopic topological density is correlated with vortices and monopoles [30]. 

The UV filtered topological density shows cluster properties reminiscent of in- 
stantons. We found a highly correlated behavior of a huge number of lowest modes 
resembling the presence of a semiclassical background. The suitably UV smoothed 
field strength tensor becomes selfdual (antiselfdual) at "hot spots" where the cor- 
responding action becomes maximal. The reason for this self-organization between 
the fermionic modes is still unclear. 
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